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We calculate the electron spectrum of semileptonic decays of B-mesons into non- 
charmed hadrons. The shape of the spectrum obtained from QCD sum rules is 
in general agreement with quark model calculations. At high electron energies the 
decrease of the spectrum is less steep than predicted by Altarelli et al. Our analysis 
yields a total decay rate T{B XuCu) = (6.8ib2.0) 10-*^^ iKftP For the spectrum 
integrated over electron energies in the interval 2.4 GeV < Eg < 2.6 GeV in the 
laboratory frame we obtain T{B Xueu) = (5.0 it 1.6) lO^'^ iKifcps"^. Estimating 
all sorts of theoretical uncertainties, we obtain \Vub\ = 0.003 it 0.001 for the b ^ u 
weak transition matrix element from the new CLEO-data. 
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I. INTRODUCTION 



Inclusive semileptonic decays of B-mesons have been studied for already more than 
a decade. Lacking a reliable identification of the hadrons in the final state, the data on 
semileptonic decays provide at present the most accurate information on the determination 
of the CKM-matrix-element \Vub\- The accuracy of the determination of \Vub\ is, however, 
limited to a great extent by the fact that the h ^ u transitions can experimentally be 
separated from the dominating b ^ c decays only in a rather narrow kinematical region of 
high electron energies, close to the kinematical threshold. 

Most theoretical predictions for the spectrum are based on a spectator model This 
approach makes use of the fact that in the limit of infinite b-quark mass the decay of 
the meson is reduced to the /9-decay of a free b-quark, which at large electron energies is 
modified by Sudakov-type radiative corrections. In practice, however, the b-quark is not 
heavy enough, and both the interaction with the spectator quark and hadronization effects 
cannot be neglected even for the total semileptonic decay rate. The numerical importance 
of l/rrih corrections becomes apparent already from simple phase space considerations: since 
the total decay rate is proportional to m%, the replacement of the hadron by the quark 
mass reduces the decay rate by nearly a factor two. Even for the ratio of semileptonic 
to hadronic width, where the dimensionful factor cancels, a calculation of the 6 — c 
transition using free quarks yields a value which is (10-15)% below the data The 
situation gets considerably worse for the electron spectrum near the charm threshold, where 
the heavy quark expansion breaks down and there is no reason to believe that Sudakov-type 
effects provide the most important correction. In existing calculations of the spectrum in 
this region the deviations from a simple free quark decay picture are taken into account 
in a model-dependent way, introducing constituent quark masses and Fermi-motion of the 
quarks inside the hadrons by a certain smearing of the momentum ||^ or by using the 
Peterson fragmentation function 0. The results show that the influence of the spectator 
quark and hadronization effects become indeed important at large electron energies, and 
quantitative predictions depend strongly on the model assumptions and on the parameter 
values. Especially unpleasant is that the kinematical endpoint is not well defined, since it 
depends on the constituent quark masses used in the model and not on the mass of the 
hadrons. 

In an alternative approach advocated in Refs. the electron spectrum of the inclusive 

decay is obtained by an explicit summation of contributions of various exclusive channels, 
e.g. B —>■ vrez/, B peu, which are calculated in the framework of a quark model. In 
this method, the results rely heavily on the assumed shape of the exclusive form factors. In 
addition, it is not possible to take into account exclusive decays with two or more nonresonant 
mesons (e.g. two pions in an S-state), which, in contrast to the decays to charmed final 
states, are expected to play a nonnegligible role [0. Recently, a two-component model 



has been suggested in Ref. [jTO[, which aims to combine the positive and avoid the negative 
features of both approaches. However, the justification of this two-component model has 
been questioned in Ref. |Tl|. 

In view of conflicting results in the current literature, we find it important to address the 
problem from a completely different point of view. In this paper we calculate the electron 
spectrum of the inclusive semileptonic B-decay within the framework of QCD sum rules. 
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This method is due to Shifman, Vainshtein and Zakharov [|T2|. It is essentially a matching 
procedure between the operator product expansion for a suitable correlation function in 
the not so deep Euclidean region and the representation of the same correlation function 
as dispersion integral in terms of contributing hadronic states. Nonperturbative effects are 
taken into account by nonvanishing vacuum expectation values of gauge invariant operators, 
the so-called condensates. QCD sum rules have proved to be a very successful tool for 
calculations of both static and decay properties of light and heavy hadrons. Although they 
often provide less detailed results than quark model calculations, they have the advantage 
of being closer to the first principles of QCD. In recent papers we have applied this method 
to the study of exclusive semileptonic decays . The application to inclusive decays 

poses some new problems, which are discussed below. We show that it is possible to use 
the sum rule technique in order to calculate the electron spectrum up to an electron energy 
Ee — 2.1 GeV in the rest frame of the B- meson, which is slightly below the charm threshold. 
The shape of the spectrum and a constraint at the endpoint, however, allow an extrapolation 
to the experimentally accessible region E^, ^ 2.3 GeV without introducing a large numerical 
uncertainty. Predictions for the shape of the spectrum and integrated rates are given and 
compared to other model calculations. 

Our paper is organized as follows. In Sec. II we recollect necessary kinematical formulae. 
In Sec. Ill we derive sum rules for the electron spectrum and discuss their region of appli- 
cability, which is constrained by the requirements of a sufficiently large interval of duality 
and the existence of the short-distance expansion. In Sec. IV we discuss the heavy quark 
limit of our sum rules. Sec. V contains numerical results and conclusions. Technical details 
of the calculation are given in the appendix. 



II. KINEMATICS 

The total width F of the inclusive decay B Xue~iJ, where Xu denotes a charmless 
hadronic state, can be written in terms of a leptonic tensor L^^, and a hadronic tensor W^,^ 
as 

^ ~ 327r^msJ2E, 2E, ^ ^^'^^ 

with 

L^^ = 2 {p^^p: + p^^p: - g^^p^p, + ie^%y,pi) (2.2) 

and 

(27r)3 5\pB -Pe-p.-J2P^)iB\ jI{0) \n){n\ j^(0) \B). 

i 

(2.3) 

The summation includes all possible n— particle final hadronic states; = ^7^(1 — 75)6 is 
the weak current mediating the decay b u, and pe, Pu, Pb are the four-momenta of the 
electron, antineutrino and B-meson, respectively. 



E/n 



=1 . 



d^Pi 
{27r)^2Ei 
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The hadronic tensor W^,y can be expressed in terms of invariant functions: 

W,. = -f^e^^^{B\[j,{x),jm]\B) 
J Zn L J 

= ^Im t Jd'x e^"^ ( B I Tjt (O)j^(x) | B ) 

= ~ g^,^Wl+pB^PBuW2-ief,^p^p''Bq''W3 + ... (2.4) 

The missing terms proportional to different Lorentz structures do not contribute to the 
decay rate in the hmit of vanishing lepton mass. The structure functions Wi depend on the 
invariant mass squared of the hadrons in the final state, s = {pb — Pe — PuY, and on the 
invariant mass squared of the lepton pair, = (pe 

Integrating over the neutrino momentum, we obtain from 



327r3 



m 



JdEe dq^ ds WWi + ml (^[^1 + - s] - 2El - M W2 



(!^£±^ _ ,E.) (2.5) 

with the electron energy Ee in the rest-frame of the B-meson. The limits of integration are 
given by: 

< i^^e < (2.6a) 

< g2 < 2mBEe, (2.6b) 

< s < (m^ - 2E,){mB - qV{2E,)). (2.6c) 

The subprocess with the lowest possible value of s is the production of a single pion in the 
final state. Hence the structure functions Wi vanish for s less than and the accessible 
phase space is restricted to the values depicted in Fig. |^. The shaded area indicates the 
region of phase space where s < 1 GeV^, the dashed-dotted hues bound regions with s < 
5, 10, 15,20GeV^ (from top to bottom). Since in the present paper the calculations are 
carried out in the chiral limit — 0, we will in the following stick to the limits of integration 
given in ( |2.6aD . 



III. DERIVATION OF THE SUM RULE 

Our objective is to evaluate the structure functions Wi entering (p.5| ) by means of QCD 
sum rules. To this end we consider the four-point function 

T,, = t'[d'xdSjd^ze-''^^^-'P^y+^P^^{0\ TjBiz)jliO)j,ix)jUy) I 0) 



Tggf,u + Tp^iPi + P2)/.(pi + P2)u + T,]^t^yp„{pi + p2Yql + ... (3.1) 
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where Jb = qi'Jbb is the interpolating field of the B- meson and = m7^(1 — 75) & is the weak 
current mediating the 6 — > m decay. 

The four-point function in can be calculated in perturbation theory in the deep 

Euclidean region, i.e. for large negative values of the virtualities s, {pi — ml) and — ni^). 
When the physical region is approached, nonperturbative effects become important. The 
idea underlying the QCD sum rules method is that at virtualities of ~ — 1 GeV^ the most 
important nonperturbative corrections can be taken into account by nonvanishing vacuum 
expectation values of operators occurring in the operator product expansion of nonlocal 
matrix elements like the correlation function ( ^.11 ). 

In our calculation, we take into account the leading perturbative contribution, which is 
given by the box-diagram shown in Fig. ^(a), and the contributions of the quark and the 
mixed condensate, examples of which are given in Fig. 0(b). 

The correlation function (|3TT| ) at t = {pi — P2Y = and := qf = {p2 — Pi — qiY 
is related to hadronic matrix elements via dispersion relations. We express the invariant 
functions Ti{s , p\, P2, q^) by the dispersion relation (see e.g. ||T6| ) 



ns.plpW)= fd~s'^^lMl]^+ L^^ U~s.pIpW) _ (32) 
J s — s J s — s 

ml —00 

The dispersive part for s > is due to the insertion of intermediate states in the di- 
rect channel, ( | JbJI I ) ( I J^i j1? I ), whereas the dispersive part below Su comes from 
the crossed channel ( | j/x I ) ( I I ) • The lowest lying intermediate state Zu 
contains two heavy quarks und thus 

Su = pI+ pl + '^q^ - {rriB + mB*)'^ . (3.3) 

In the heavy quark limit rriQ — * 00, q'^/rriQ < 1, the contribution of the crossed channel 
vanishes, and thus the structure functions can be related directly to the expansion at short 
distances @. On the other hand, for finite quark masses the two cuts may come very 
close to each other for on-mass-shell B- mesons and — m^. It is specific to the QCD 
sum rule approach, however, that calculations are done for sufficiently off-shell values of p^: 
pf — ^ — 1 GeV^. In this way the singularities of the direct and the crossed channel are 
well separated. Furthermore, the application of this technique requires to stay several 
GeV^ below m^, the quark mass squared, and hence the branching point lies at least 
6 GeV^ below the the physical region of the B-decay. We shall see later that in our approach 
the contribution of the crossed channel can be neglected (see below Eq. ( |3.8| )). The required 
information on the structure functions Wii^s^q^) of the B-decay is contained in the terms 
with poles in both the variables p\ and p^ of the functions Wi{s^p\^p\^ q^) in Eq. ( |3.2|) : 



Here the dots stand for contributions of higher resonances and the continuum. 

In order to estimate the contribution of higher states in ( |3.4D we use an approximation 
which is standard in the QCD sum rule approach and assume that all contributions of 
higher mass states are eliminated by expressing the four-point correlation function in (p.l|) 
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by a double dispersion relation in and and retaining only the contribution below a 
certain threshold Sq. The interval of duality, i.e. the value of the continuum threshold Sq, is 
estimated by the criterium of stability of the sum rules. 

Following [|l^], we subject the correlation function to a Borel transformation in 

both the variables pi and P2- For an arbitrary function of the Euclidean momentum, f{P'^) 
with = — p^, this transformation is defined as 

p2/N=A'P fixed 

where is a new variable, called Borel parameter. The objective of this step is to reduce the 
dependence of the sum rules on the unknown higher-order terms in the operator expansion 
and on the continuum model. Indeed, the application of the Borel transformation to a 
typical contribution yields 

(p2_^2)n ^ (M2)- ' ^ ' 

and thus contributions of vacuum condensates of high dimension (which contain high powers 
of (p^ — m^) in the denominator) are suppressed by factorials. In addition, the contribu- 
tions of higher mass states become exponentially suppressed, since the Borel transformed 
expressions for the hadron propagators (p^ — m^)~^ and (p^ — m^)"^ differ by an exponential 
factor exp(— — mD/M"^). It is important to note that a single Borel transformation in 
the variable pf = pi = is not sufficient for our purposes, since it does not suppress terms 
containing a pole in one variable only - either in pi or in P2. 

Since, according to Eq. ( p.2| ), we are further interested in comparing T^,^ to the dispersion 
integral in s in terms of the structure functions Wi and Ri, we express the amplitudes Tj in 
form of triple dispersion integrals in s, p\ and p^. This step requires some care, since in the 
region of positive the condition t = (pi — ^2)^ = enforces p\ = p^. Thus, we start from 
negative values of g^, where pf and pi are independent variables and write 

Ti{plpls,q^,t = 0) = [d~sdsids2- ^iw'''''!'2wg + 

As can be shown, ( |3.7| ) defines an analytic function of and can be continued to positive 
values of q^. As pointed out in [|13|, for positive q^ non-Landau singularities may generally 
enter the game. These additional contributions are, however, numerically unimportant in 
the present analysis and thus have been omitted, cf. the appendix for details. 



Due to the condition t = 0, the support of the spectral function pi in (|37^ ) is restricted 
to the line si = S2, i.e. pj(si, S2, 5, g^) = pi(si, s, — S2), which means that the four- 

point function does not receive any contribution from non-diagonal transitions like B B' . 
All contributions to the four-point function with a single pole only in p\ or pi cannot be 
expressed as double dispersion relation in both p\ and pi and are related to subtraction 
terms in one of the variables. All these subtractions are eliminated by the double Borel 
transformation. 



Putting together (|3.1|), (p.2|) and (|3.7|) we get, after performing the Borel transformation: 
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K0b.|5)Pe--i/-^J fd~s^^^^+ fd§?^\= fdsjds/^^''^'^'^' 

li/ S S "J S S \ "J "J 

\ml -oo Jo 

(3.8) 

where Ri{s,q^) is defined analogously to Wi{s,q^) in (^.41). 

Staying strictly in the framework of the operator product expansion we cannot discrim- 
inate distributions from the direct and the crossed channel. We, therefore, have to invoke 
a concept of duality and assume that the relation ( p.8| ) holds true for the spectral densities 
themselves, if the latter are smeared over a sufficiently large interval of s. The leading con- 
tributions to the operator product expansion of ( p.7| ) do not contribute to the discontinuities 
Ri of the crossed channel; for this reason the crossed channel singularities do not occur in 
our calculations. 



Existing calculations of deep inelastic structure functions by QCD sum rules p!8|Jig 
use the approximation of local duality, equating in the expressions analogous to ( |3.8D the 
structure functions and the spectral densities pointwise in s. Since we are interested in the 
electron spectrum and integrate over the neutrino momentum, we automatically are led to 
an integration over s, cf. Eq. ( |2.5| ). Therefore, in this paper we only make the much weaker 
assumption, that the physical spectral density coincides to the one obtained by the operator 
product expansion upon the integration over a sufficiently large interval of duality As as 
mentioned before. Since the lowest mass physical state is the pion, we demand Smax to be 
larger than the interval of duality in the QCD sum rules for the pion As ~ 1 GeV^ [112 



Since at Ee, k, Emax the physical region in s shrinks to a point, the requirement of having 
a sufficient interval of duality excludes a certain fraction of the kinematical region, which 
we show as a gray area in the Dalitz-plot Fig. |l|. For small values of E^. it covers only a 
small fraction of the integration region in and is heavily suppressed by phase space, so 
that this limitation is not restrictive. More important constraints on the accessible values 
of follow for higher electron energies. However, the major part of the excluded region in 
the Dalitz plot is not accessible to the sum rules for a different reason as well: the value of 
q^ — ml should be sufficiently large and negative in order the operator product expansion 
to be justified. Taken altogether, these requirements lead to an upper bound for accessible 
electron energies of 2.1 GeV. At this energy requiring an interval of duality larger than 
As ~ 1 GeV excludes 4.7% of the whole phase space in s and g^. At lower electron energies 
it is an even smaller fraction. 

Expanding (|3.8|) in powers of energy one may relate the moments of structure functions, 
/ dss~^Wi{s), to corresponding integrals over the spectral densities, and, in turn, to matrix 
elements of local operators over the B-meson, which appear in the short-distance expansion 
of the T-product of weak currents. The operator expansion of this product is the starting 
point in approaches to the inclusive B-decays based on the l/rrib expansion or the parton 
model |T^. It is worthwile to note that the integration over s in the differential decay rate 



at fixed Ef. in ( p.5|) is constrained by the value 

Smax = - 2^e) (^B " ^ j (3.9) 

which for not too large E^ is less than the actual support of the structure functions. Thus, 
the knowledge of the first few moments of the structure functions Wi is not sufficient for the 
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calculation of the inclusive spectrum, and we cannot reduce our task to the calculation of 
three-point instead of four-point functions. 

Writing the lepton decay constant ( | I -S ) as m'^sfB/^b, we get from ( p.8| ) the fol- 
lowing sum rules for the structure functions Wf 



9 *0 

^ii^^Q^) = 4 % dsipi{si,q^,s)e 
^bJb 



-2(si-m|)/M2 



(3.10) 



These equations should be understood in the sense that the physical Wi coincide with the 
ones given in ( p.lO| ) after smearing over an interval of duality in s larger or of the order of 
1 GeV2. 

The spectral densities pi receive contributions from perturbation theory and the vacuum 
condensates 



p^ = p^i^''^ + pf {qq)+ pf ( qcrgGq). 
A calculation of the box graph in Fig. |^(a) yields (see appendix for the details) 



pi 
Pi 



'crt 



3 ml — si 

4v]-2 g2 



m. + 



2q^s + {ml - g^)(s + - SiY 



9{si - s 



^pert 



pPert 



3 ml — Si 
27r2 AV2 



where 



T 2 smt 
s[ = ml+ , " 
m^ — g2 

Note that pf^'^* is regular in the limit — 0. 

The contribution of the quark condensate is given by 



(3> 
Pi 



(3> 
P2 



= -2mi,{s + ml-q^)5{si-ml)5{s), 
= —8mb 5{si — ml) S{s), 
pf^ = 4mb5(si - ml)6{s), 
and the contribution of the mixed condensate reads 

5{s,-ml)-\2iml-q^)5is) 



(5) 



-6is) + (ml - q')5'is) 



liml-qr5'is) 



[si - H b[s)d [si - 



(3.11) 

(3.12a) 

(3.12b) 
(3.12c) 

(3.13) 



(3.14a) 
(3.14b) 
(3.14c) 



(3.15a) 
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) + \-A6{s) + 2{ml - q^)6\s)\ 6'{s^ - ml) 



+ 2ml5{s) 5'\s,-ml)} 



rrih \^S'{s)S{si 



(3.15b) 



(3.15c) 



Note that 6{s) has to be understood as the hmit of 6{s — mD where we have put m^, the 
mass of the u-quark, to zero. Thus there is no ambiguity in the integration over s within 
the range specified by ( p^.bal ). Contributions of the u-quark condensate (mm) and of the 
corresponding mixed condensate (ucrgGu) vanish. Note that the contributions containing 
6'{s) give rise to terms ~ 6{2Ee{mB — m^/ {ms — 2Ee)) — q^) in the double differential rate 
(fV/{dEJq^). 

The leptonic decay constant Jb entering Eq. (|3.1CI| ) is determined by a two-point sum 
rule where the same approximations are made and the same input parameters are used 
as for the sum rule (p.lO|) , except for the value of the Borel parameter M. 

MV2 |T3|: 



two times smaller, M|pj = 



2pf, which is taken 



(3.16) 



3ml 
8^ 



so 



dsi 



m 



2\2 



si)/M? 



2pt 



m, 



mt 



2 



2pt 



(qagGq). 



Considering the ratio of sum rules instead of fixing Jb at a certain value considerably reduces 
the dependence of the results on the input value of the b-quark mass and on the value of 
the Borel parameter, since the effect of a change in these parameters is cancelled between 
the numerator and the denominator to a great extent. 

Radiative corrections to the sum rule in ( 3.10 ) may in general be important, cf. ||2l|,22]. In 
the limit of large m^, however, all radiative corrections factorize into corrections to Jb, which 
are cancelled by taking the ratio of the four-point and two-point sum rule, and Sudakov- 
type radiative corrections to the spectrum of the decay of a free b-quark, which we take into 
account as a multiplicative factor, see next section. A full account for the radiative correction 
requires a laborous calculation of the a^-correction to the box graph and is beyond the tasks 
of this paper. On formal grounds, all corrections which we do not take into account are 
suppressed by a power of the heavy quark mass. 



IV. THE HEAVY QUARK LIMIT 



We consider the region where the invariant mass of the lepton pair and the electron 
energy 2Ee constitute a finite fraction of the heavy quark mass, i.e. 



m^-q 
mb - 2Ee 



I), 



0{m_ 
0{mt). 



(4.1) 
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It is convenient to introduce the scaling variable 0] 



X = 2Ee/mB (4.2) 

in terms of which the last of the conditions in (|4.1| ) reads 1 — x = (9(1). In this kinematical 
range the sum rules ( |3.10| ) simplify drastically in the limit of infinite heavy quark mass and 
yield the results of the free quark decay. 

To show this, we use the common assumption that all excitation energies remain finite 
in the heavy quark limit, and hence sq, the continuum threshold, scales as 

So-m2 = (mB + ^exc)^-m2~0(mf,). (4.3) 

The working region in the Borel parameter is determined by the requirement that the ex- 
ponential suppression factor for the contributions of excited states, exp[— (sq — m%)/M^], 
remains finite in that limit. Therefore, also the Borel parameter should be taken at 
values of order 0{mb). Next, we note that the support of the structure functions is re- 
stricted to an interval < s < 0{mb). As far as the perturbative contribution to the sum 
rules is concerned, this statement follows from the kinematical restriction for the double 
discontinuity of the box graph, Si > Sl, in ( |3.13|) , where from 

s<{s,-ml)^^:^^^<0{m,). (4.4) 
mi 



The last inequality is a consequence of Eq. (|4.3|) . Therefore s can be neglected against si, 
which is of 0{ml). The nonperturbative contributions generally form a series in (5-functions 
of s and its derivatives, which should in principle produce a certain smooth function upon 
summation. 

Retaining the leading terms in the heavy-quark mass only, we find from (|3.10|) : 



^W,{s, q-) = ds, e^(-|--)/-^^(., - ml 

- 2{qq)mb{ml - q^)e^^<-<'^"''" 5{s) 



ml 



+ 2{qagGq)^^{ml - q')e'^<-<y''^' 6(3) (4.5) 

In order to obtain the inclusive electron spectrum, we integrate ([4.5| ) over s. Changing the 
order of integration in s and in si, we find the integration intervals 

ml<si< So, 

0<s<mtn (si - ml), {mi, - 2Ee)(mb - ^) ) . (4.6) 



Due to ( [4.4| ), the second term in min(. ..,...) is parametrically large as compared to the 
first one. Hence the full support of Wi in s is covered, exemplifying the general statement 
that in the heavy quark limit the knowledge of the first moment / dsWi{s) is sufficient for 
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the calculation of the electron spectrum in the kinematical range given in ( |4.1|) , cf. Refs. 
Thus, we obtain: 



ml 



{3 f^o 



2\2 



ml 



qq )mbe 



2{m%-ml)/R4^ 



m 



+ {qcT9Gq)j^e 



b2iml-ml)/R4^ 



(4.7) 



Letting = 2M^^^, one realizes that the expression in curly brackets coincides (apart from 
an overall factor ml) with the leading term of the two-point sum rule ( p.l6| ) in the limit of 
infinite quark mass. The same factorization takes place in the other two structure functions. 
Thus, in the heavy quark limit we end up with the following expressions: 



dsWi{s,q^ 



2{ml - q') 



J dsW2{s,q^) = 8, 



(4.8a) 



(4.8b) 



dsW3{s,q^ 



(4.8c) 



which coincide with the structure functions for the decay of a free heavy quark. The corre- 
sponding electron spectrum is p| 



dT 

dx 



G^\Vut\'ml 
967r3 



x^{3-2x), 



where x = 2Ee/mi,, and the total decay rate equals 

^ _ G'\V^,\^ml 
1927r3 



(4.9) 



(4.10) 



Note that to this accuracy there is no distinction between the quark and the meson mass. 

The results presented above do not take into account radiative corrections oc as{mb). As 
it is well known, such corrections can in general get enhanced by large logarithms of the 
heavy quark mass. The corresponding corrections to the spectrum for the decay of a free 
quark have been calculated in H] . They result in a modification of the uncorrected spectrum 



(4.9) by a multiplicative factor: 



dT 

dx 



dx 



2a, 
37r 



G{x) 



exp 



-^ln^(l 
37r ^ 



X] 



(4.11) 



where the Sudakov-like exponential has been factored out and the function G is given by 



G{x) 



X 



ln(l — x) 



3, l + 2x 11 

- m 1 

2 3 6a; 



Ax 



2 



:i-x) 



2 V3 



X [ —Tx'^ 

\3 12 



+ 2x Li {x) 
19\ 11 



6 



(4.12) 
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In numerical calculations, the results of which are presented below, we take into account 
this Sudakov-type correction, which is the leading one in the kinematical region in ( |4.1D . 
All the radiative corrections due to the spectator quark are suppressed by a power of m^. 
It should be noted, however, that for the region of high electron energy 1 — x ~ l/m^ the 
application of ( |4.11| ) is not justified, and the true effect of radiative corrections may be 
different, as indicated by the qualitative change of the Sudakov-type exponential behaviour 
with the account of the constituent mass of the u-quark [H]. 



V. RESULTS AND DISCUSSION 

The results of the numerical calculations given below have been obtained using the 
following values of the vacuum condensates: 

(gg)(lGeV) = (-0.24GeV)^ 

{qagGq){lGeV) = 0.8 GeV^ ( GeV). (5.1) 

These values have been rescaled to a higher normalization point = M^, which is the 
relevant short distance expansion parameter in the sum rule (|3.10|) . We have varied the value 
of the pole mass of the b-quark rrif, within the range (4.6 - 4.8) GeV, which covers the values 
given in literature |^,24|. The continuum threshold Sq is determined from the requirement 



of stability of the two-point sum rule in ( |3.16| ) and depends on the quark mass. We use the 
values soimt, = 4.6 GeV) = 36 GeV^ So(4.7 GeV) = 35 GeV^ and So(4.8 GeV) = 34 GeV^ 
7^. These values provide a good stability of the two-point sum rule in a large range of 



Borel parameters M|p^. Note that we use the same value of sq in both the four-point and the 
two-point sum rules. Similarly, the range of appropriate values of the Borel parameter 
in the four-point sum rule is fixed from the two-point sum rule ( p.l6| ) by the requirement 
that both the contribution of the continuum and of nonperturbative corrections do not 
exceed 40% each. In this way we come to the range 3.5 GeV^ < Mlp, < 5 GeVl Similar 



to the case of the three-point function [|13|, the Borel parameter of the four-point function, 
M^, should be taken two times larger than M|p^, i.e. we arrive at the interval of values 
7 GeV^ < Af2 < 10 GeVl 

The QCD sum rule predictions for the differential spectrum dV/dx in the rest frame of 
the B-meson are presented in Fig. ^ in dependence on the scaling variable x = 2Ef./mB for 
different values of the Borel parameter, and for input values mb = 4.8 GeV and Sq = 34 GeV. 
The curves are calculated by means of Eq. ( p.5|) up to the electron energy 2.1 GeV (x = 0.80), 
where the upper limit of integration in is still small enough as to ensure the validity 
of the operator product expansion (cf. Sect. III). For higher electron energies we did a 
smooth interpolation to the end point of the spectrum at x = 1, where dT/dx = owing to 
vanishing phase space. The position of the maximum of the spectrum is obtained at values 
X ~ 0.75, which are still within the region of validity of the QCD sum rule calculation, 
and due to that the uncertainty introduced by the degree of the interpolating polynomials 
is negligible: it turns out to be smaller than 1% in the experimentally interesting region 
2.2 GeV < Ef. < 2.4 GeV when changing from a cubic to an order five polynomial. We 
find that the variation of the Borel parameter within the range given above induces a ±15% 
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uncertainty in the absolute normalization, while the shape of the spectrum is much more 
stable. 

Fig. ^ shows the separate contributions to the differential spectrum coming from the per- 
turbative graph and from the nonperturbative corrections. The contributions of the vacuum 
condensates are large by themselves, but partly cancel in the sum. Up to the very end of the 
spectrum the relative weights of the perturbative and nonperturbative contributions remain 
practically constant. 

A variation of the mixed condensate in the range 0.6 GeV^ < ttiq < 1 GeV^ leads to a 
change of the decay rate by ±10%. 

The influence of the mass of the heavy quark, which is large in the case of the sum rule for 



/b [PIIkSI , is eliminated to a great extent by compensations between the numerator and the 
denominator in ( p.lO| ). A variation of from 4.8 to 4.6 GeV reduces the decay rate by 8%. 



A similar compensation occurs for the dependence on the continuum threshold Sq: variation 
between 30 and 38 GeV^ results in a change by ±4%. Variations of other parameters within 
the standard limits have only marginal effect. 

The radiative corrections in ( |4.11| ) reduce the decay rate by 20%. Since the main effect 



is however achieved at large values of x, where neither the derivation of ( [4.11| ) nor the QCD 



sum rules approach are justified, we add an additional error of 10% to our results, which 
corresponds to an uncertainty in the radiative corrections of 50%. 

It is worthwhile to mention that major part of uncertainties of the QCD sum rule calcu- 
lation result in an uncertainty in the overall normalization factor, but not affect the shape 
of the spectrum. 

Combining all sorts of errors, we obtain the total decay rate: 

r = (4.5 ± 1.3) lO'^^lK^r GeV, 

= (6.8±2.0) 10^^ iKfeps^^ (5.2) 

In order to be able to compare to the experimental results which were obtained for B-mesons 
created on the T(4S') resonance, we have to boost to the laboratory frame, 

where the electron energy is measured in the laboratory frame and vb = 0.065 is the 
velocity of the B-meson. The change induced by the boost is negligible for the total rate, 
but modifies the high energy part of the spectrum. Thus, in Table I we give the values of 
the integrated spectrum in energy bins of 0.1 GeV from 2.1 GeV to 2.7 GeV both in the 
rest and the lab-frame. 

For the integrated value of the decay rate in the interval of electron energies 2.4 GeV < 
Ee < 2.6 GeV we get: 

r''** = (5.0 ± 1.6) 10^2 \Vub\^s'\ (5.4) 

Combining the latter value with the new experimental result of CLEO [0] for the branching 
ratio in the same region of energies. 
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B= (0.53 ± 0.14 ± 0.13) 10-^ 



(5.5) 



we obtain the following value for the CKM- matrix element: 

iKbl = (0.003 ±0.001) (^^^^ ^ , (5.6) 

where tb is the lifetime of the B-meson. The error takes into account the accuracy in the 
calculated decay rate, a conservative estimate for the uncertainty induced by the extrapo- 
lation and the experimental error, where theoretical and experimental error are roughly of 
the same size. 

The comparison with the other existing calculations of absolute rates is given in Table I, 
and the comparison of the predicted shape of the spectrum is shown in Fig. ^. The QCD 
sum rule calculation of the electron spectrum for 2.1 GeV < -Eg"* < 2.5 GeV lies inbetween 
the predictions of the free and the modified P,p!0| and the direct summation of the 
contributions of exclusive channels 0, while for higher energies our results coincide with 
and |T0] within the errors. At such high energies, however, the extrapolation procedure used 



might not be adequate. All calculations of the shape of the spectrum are in a good agreement 
at low values of x, while starting from x ~ 0.5 the discrepancies amount (10-15)%. The 
calculations using a certain modification of the free quark model have an endpoint which 
might differ from the physical one, since the quark masses enter in this approach instead of 
the hadron ones. The slope of all curves near the end is, however, not so much different, the 
one of being the smallest. 

Concluding we remark that the sum rule method allows a determination of the semilep- 
tonic decay rate which is on a theoretically sound basis for electron energies below 2.1 GeV. 
Though at first sight it may seem unsatisfactory that the most interesting range of electron 
energies can be reached only by a constrained extrapolation, this does not induce a major 
error. This is due to the fact, that the electron spectrum can be calculated to energies larger 
than the position of the maximum and that the endpoint of the spectrum, where it vanishes, 
is determined by hadron, and not quark masses in our approach. In the determination of the 
CKM matrix element Vbu the uncertainty of the lifetime is still the most significant error. 
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THE TRIPLE DISPERSION RELATION OF THE BOX-DIAGRAM 

In order to be able to subtract the contribution of the continuum from our sum rules 
and to take the imaginary part in s, the invariant mass squared of the final state hadrons, we 
need to represent the box-diagram Fig. ^(a) by a triple dispersion relation. Since the forward 
scattering amplitude is relevant to our problem, we consider the kinematical configuration 
t = {pi — = {qi — = 0. In addition to that, we want to keep pi ^ p\ whilst 
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qf = q2 = which imphes q^ < 0. Thus, in our case the independent variables determining 
the box-diagram are p\, q^ and s. The external momenta of the interpolating field of the 
B-meson, pi and p2-, are considered in the Euclidean region p\.,p\ < 0. 

We infer the existence of such a triple dispersion relation for the fermionic case from the 
analysis of the scalar box-diagram for < 



S 



(27r)4 ((pi + kf - m2)((p2 + kf - ml){p, + q^fk"^ 

5{si - S2) 



Ids- (dsi [ds 

J s — s J J 



IGtc'^ ml — q'^ J s — sJ J {si — Pi){si — P2^ 



IBtt^ ml — q^ J^^ s — s J "^^ (si — pl){si — pi) ^'^'^ 

sf 



with 



^f = 4^+< (8) 
m^ — q'^ 

The right hand side of (|^ can be continued analytically to > with p\ and p\ remaining 
independent variables. We have checked this issue also explicitly by comparing with the 
corresponding Feynman-parameter integral. For arbitrary values of g^, Eq. (0) can equally 
well be obtained from the scalar triangle graph T with external momenta squared p\, q^ and 
s as 



s = ^npis,q^) 



(9) 

(Sl-p?)2^{si-pf){si-p2) 



which follows immediately from the expressions for the loop-integrals. 

As pointed out in |]l3l, a naive continuation of the double dispersion relation for the 



triangle graph to the region of positive values of q^ is not possible because of the existence 
of non-Landau singularities. Technically, the problem arises because the double spectral 
function of the triangle graph T with respect to p\ and s, pT{Pi,q'^, s), contains a factor 
1 / A^/^, where X = pf + s"^ + q^ — 2pls — 2p^g^ — 2sg^. For positive values of g^, A can vanish 
and pt can become singular at the boundary of the integration region as given by the Landau 
equations. This allows the logarithmic branching cut on the unphysical sheet of the square 
root to dive through the square root cut into the physical sheet. The proper treatment 
of that additional singularity is described in detail in [|13| . As follows from Eq. (|^, for 



the scalar box-graph the non-Landau singularities are absent, but are present for fermions. 
They play, however, no role in our analysis, because it turns out that for physically relevant 
g^ they lie above the continuum threshold or are numerically negligible due to phase space 
suppression. 

It proves convenient to apply the above method of taking the derivative of some triangle 
diagram with appropriate vertices with respect to ml likewise to the fermionic case with 
open Lorentz-indices. From the double spectral functions for the tensor-integrals given in 
11 we obtain (for g^ < 0) 
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(27r)4 ((pi + kf - ml){{p2 + kf - ml){pi + qi)^k^ 



P+^ipi + P2)^^ + P-^ipi -P2)M + QigiM- (10) 



The only invariant needed is P+ for which we obtain: 



oo 







with \ = sf + s"^ + q'^ — 2sis — 2sig^ — 2sg^. 

The integral with two open Lorentz-indices is decomposed as 



(27r)4 ((pi + kf - ml){{p2 + kf - ml){pi + qi^k^ 



Ag^, + B^{pi+P2)^{pi+P2). + ... (12) 



with 



t 1 {ml - g^)A^/^ + 2g^g + (mg - g^)(g + q^ - Si) 

IQT^^Aq^n-s-sj'^'' X^I\s,-pI){s,-pI) ^^^^^ 



oo ^ oo 

B — Ids- Idsi 



IGTT^y s-sJ {si - pI){si - pI) 
r 1 , 2 , K-g^)(g^ + g-gi)+2gg^ | 

which are the only invariants needed. 

The expressions given above (including non-Landau singularities for q^ > 0) have been 
checked for arbitrary values of q^ by comparison with the corresponding Feynman-parameter 
integrals and by an independent calculation of the double spectral densities using the de- 
composition of momenta in light-cone variables. 

By means of the above expressions we obtain the perturbative double spectral functions 
given in Sec. III. 
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FIGURES 



FIG. 1. Dalitz-plot of the phase space in the plane of invariant lepton mass squared, g^, and 
electron energy, Ef>. The solid border-line where s = can only be reached with a single pion in 
the final hadronic state. The shaded area shows the region of phase space with s < 1 GeV^. The 
dashed-dotted lines denote curves of constant s > 5, 10, 15, 20 GeV^ (from top to bottom). 



FIG. 2. Feynman diagrams contributing to the operator product expansion of ( |3.lD : (a) per- 
turbative contribution, (b) contribution of the quark and part of the contribution of the mixed 
condensate. Lines ending in a cross denote vacuum expectation values. 



FIG. 3. The spectrum iVubl^"^ dT / dx as function of x with mf, = 4.8 GeV and sq = 34 GeV^ for 
different values of the Borel parameter, = 7,8,9,10 GeV^ (from top to bottom). Solid lines: 
evaluation using ( 3.1C| ), dashed lines: interpolation to zero at j; = 1. 



FIG. 4. Relative contributions of the separate terms in the operator product expansion to the 
spectrum as function of x (M^ = 8 GeV'^, nib = 4.8 GeV, sq = 34 GeV^). Solid line: perturbative 
contribution, long dashes: contribution of the quark condensate, short dashes: contribution of the 
mixed condensate. 



FIG. 5. The spectrum \Vub\~'^dT /dx as function of x with mi, = 4.8 GeV, sq = 34 GeV^ and 
= 8 GeV^. Solid line: the spectrum using ( 3.10| ), long dashes: the same with inclusion of 
radiative corrections for the free-quark decay acc. to ( [4.11| ). 



FIG. 6. 1/TdT/dx as function of x for the free quark decay and in different model-calculations. 
The chosen normalization emphasizes the differences in shape. ACCM: GISW: Q, RDB: |1C|, 
BBD: this paper. 
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TABLES 



TABLE L Central values for the inclusive semileptonic B e v decay rate from this 

paper and from different other calculations. Free: free quark decay with QCD corrections 
Eq. (25), as = 0.24, rrib = 5.0 GeV; ACCM: decay including bound-state corrections [|], Eq. (27), 
nis = O.lGGeF, pf = Q.lbGeV,as = 0.24; GISW: [|], Fig. 6; RDB: [|lO|, Fig. 2, cut-off 1.5 GeV; 



BED: this paper. We give values for the total rate F in |l^f,plO^^ s~^ and for the spectrum inte- 
grated over electron energies in bins of 0.1 GeV in the rest and in the lab frame. In the lab frame 
the B mesons are assumed to have a momentum of 0.34 GeV and isotropic angular distribution. 



Frame 


Model 


F 


2.1-2.2 


2.2-2.3 


2.3-2.4 


2.4-2.5 


2.5-2.6 


2.6-2.7 


Rest 


BBD 


0.68 


0.053 


0.050 


0.044 


0.034 


0.019 


0.002 




Free 


0.85 


0.066 


0.067 


0.062 


0.045 










ACCM 


0.89 


0.070 


0.071 


0.068 


0.052 


0.019 







GISW 




0.033 


0.030 


0.023 


0.013 


0.003 







RDB 


0.92 


0.071 


0.070 


0.066 


0.043 


0.009 





Lab 


BBD 


0.68 


0.052 


0.049 


0.042 


0.032 


0.018 


0.008 




Free 


0.85 


0.066 


0.066 


0.058 


0.037 


0.017 


0.0015 




ACCM 


0.89 


0.069 


0.069 


0.063 


0.045 


0.025 


0.008 




GISW 




0.032 


0.029 


0.022 


0.013 


0.006 


0.002 




RDB 


0.92 


0.071 


0.069 


0.059 


0.039 


0.019 


0.005 
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